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LATERAL AERODYNAMIC PARAMETERS EXTRACTED 
FROM FLIGHT DATA FOR THE F-8C AIRPLANE 
IN MANEUVERING FLIGHT 

William T. Suit 
Langley Research Center 

SUMMARY 

Flight -test data have been used to extract the lateral aerodynamic parameters of the 
F-8C airplane at moderate to high angles of attack. The data were obtained during pertur- 
bations of the aircraft from steady turns with trim normal accelerations from 1.5g to 3.0g. 
The angle-of-attack variation from trim was negligible. 

Although wind-tunnel data indicate that the rolling and yawing moments are somewhat 
nonlinear with angle of attack, the angle-of-attack variations are small; therefore, the 
linear aerodynamic coefficients extracted from the flight tests permit computation of motion 
time histories which are in close agreement with the measured time histories. The aero- 
dynamic coefficients extracted from flight data were compared with several other sets of 
coefficients, and the extracted coefficients resulted in characteristics for the Dutch roll 
mode (at the highest angles of attack) similar to those of a set of coefficients that has been 
the basis of several simulations of the F-8C. 

INTRODUCTION 

The National Aeronautics and Space Administration (NASA) is currently involved in 
research on fly-by-wire control systems for aircraft. A discussion of the NASA fly -by- 
wire program is given in reference 1. The aircraft currently used by NASA as a test bed 
to study digital fly -by-wire systems is an F-8C airplane with a standard airframe. Pre- 
viously determined aerodynamics of the subject airplane came primarily from wind-tunnel 
tests and analytical calculations, and the mathematical aerodynamic model of the airplane 
was considered to be reasonable, especially for trimmed level flight. To substantiate the 
existing linear aerodynamic model for the F-8C at moderate to high angles of attack, some 
flight tests were made with the angle of attack as high as 16°. 

A maximum-likelihood extraction procedure was used to analyze the flight data. In 
this procedure, a set of equations of motion is used to calculate aircraft response to speci- 
fied control inputs. Initial estimates of aerodynamic parameters (either from theory or 



from wind-tunnel tests) are used for the initial motion computations. An iterative digital 
computer program then selects a set of aerodynamic parameters that minimizes the differ- 
ence between the computed profiles and the flight profiles. This program has been used to 
determine the aerodynamic parameters for several aircraft in the lg trimmed -flight condi- 
tion. (See refs. 2 to 4.) The details of the program are contained in reference 5 and in the 
appendix of this report. The program has not been used previously for lateral flight data 
taken at moderate to high angles of attack with a trimmed normal acceleration greater than 
lg. The program can be used as long as angle-of -attack variations from trim are small, 
so that the assumption of linear aerodynamics will not be violated. Analytical and wind- 
tunnel studies have shown that most of the lateral aerodynamic parameters can have non- 
linear variations with angle of attack over the range used in this investigation and that some 
of these variations can be significant (ref. 6). For each of the individual flights used in this 
investigation, the variations in angle of attack from the trim angle of attack were less than 
1° during 90 percent of the time history and always less than 2°. Therefore, a linear 
model for the aerodynamics was considered adequate to describe the motion of the airplane. 
The linearity assumptions were successfully used in extraction of the longitudinal aerody- 
namic parameters for the F-8C at moderate to high angles of attack, and the results are 
given in reference 7. 

The purpose of this paper is to present the lateral aerodynamic parameters of the 
F-8C airplane as calculated from flight data obtained near trim at Mach numbers of 0.7 and 
0.8, with normal accelerations of 1.5g to 3.0g. Also presented are the equations used and 
additional information on the flight data. These are followed by some comments on the 
extraction procedure and a discussion of the results of the study. 

SYMBOLS 

Values are given in both SI and U.S. Customary Units. The measurements and cal- 
culations were made in U.S. Customary Units. 

ay acceleration measured along Y body axis, g units 

b wing span, m (ft) 

c wing mean geometric chord, m (ft) 

g acceleration due to gravity, 9.81 m/sec^ (32.2 ft/sec^) 

h altitude, m (ft) 
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I moment of inertia, kg-m^ (slug-ft^) 

M Mach number 

P period of oscillatory motion, sec 

p rate of roll, rad/sec or deg/sec 

q rate of pitch, rad/sec or deg/sec 

r rate of yaw, rad/ sec or deg/sec 

S wing area, m^ (ff2) 

t time, sec 

u component of velocity along X body axis, m/ sec (ft/sec) 

V aircraft total velocity, m/sec (ft/sec) 

v component of velocity along Y body axis, m/sec (ft/sec) 

W aircraft weight, N (lb) 

w component of velocity along Z body axis, m/sec (ft/sec) 

X,Y, Z body coordinate axes through airplane center of gravity 
a angle of attack, rad or deg 

/3 angle of sideslip, rad or deg 

6 a aileron deflection (positive for left roll), rad or deg 

S g tail-plane deflection (positive for trailing edge down), rad or deg 

5 r rudder deflection (positive for trailing edge left), rad or deg 

£ damping ratio 
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6 


pitch angle, rad or deg 


p air density, kg/ (slugs/ft^) 

<j) bank angle, rad or deg 

Coefficients and derivatives: 

C L rolling-moment coefficient 

C n yawing-moment coefficient 

Cy side -force coefficient 



4 





Subscripts : 
c computed 

f measured in flight 

t trim conditions 

X ? Y,Z body coordinate axes through aircraft center of gravity 

A dot over a symbol signifies a derivative with respect to time. 


5 



EQUATIONS OF MOTION 


The equations of motion used in this study are referred to the body -axis system 
shown in figure 1 and are as follows : 

Y -direction: 


• 1 2 g 

v = pw - ru + g cos 9 sin (p + ^pV S ^ 
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Auxiliary equations: 
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iy = — (v + ur - wp - g cos 9 sin 
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(f> = p + (q sin (p + r cos <p ) tan 9 


The three-degree-of-freedom equations (eqs. (1) to (3)) were solved for the lateral 
motions. The longitudinal variables u, w, q, and 9 used in these equations were the 
flight-measured values, and hence the nonlinear contributions of these terms were used in 
the equations. 


FLIGHT DATA 
Description of Airplane 

The airplane used, a modified prototype F-8C, has been a flight -test vehicle since its 
manufacture in 1958. The F-8C is a single-seat, high-performance airplane with a single 
jet engine embedded in the fuselage. Pitch control is achieved with a unit horizontal tail. 
The center of gravity was at 29.0 percent of the mean geometric chord. The X body axis 
was parallel to and 10.16 cm (4 in.) above water line 100. (See fig. 1.) The geometric 
characteristics of the airplane are given in table I. 

Flight Tests 

The data used in this report were obtained from flights made at the Hugh L. Dryden 
Flight Research Center, and are shown in figures 2 to 6. The pilot was instructed to fly a 
coordinated turn at nominal Mach numbers of 0.7 and 0.8 with nominal trim angles of attack 
of 9° and 13°. The worst case of angle-of-attack variation is shown in figure 3, where, as 
can be seen, the angle of attack varied from the trim angle of attack by less than 2°. A 1° 
variation in a represents a 0.2g variation in normal acceleration. The airplane stability 
augmentation systems were off during the tests. The actual test conditions for each indi- 
vidual run are given in table II. The airplane mass and moments of inertia listed in table n 
were calculated as a function of the percent of fuel in the airplane and were obtained from 
tables supplied by the Dryden Flight Research Center. The mass and moments of inertia 
used are average values for the test duration. Since the mass varied less than 3 percent 
and the moments of inertia by less than 1 percent, these variations were not accounted for 
in the parameter estimation. 

Pertinent data recorded during the flight tests included lateral acceleration ayj 
the difference between total pressure and static pressure measured on a nose boom extend- 
ing 1.83 m (6 ft) in front of the airplane; pitch attitude 9; bank angle cf>; pitch rate q; 
yaw rate r; roll rate p; indicated angle of attack a and indicated angle of sideslip /3 
measured by vanes on the nose boom; pressure altitude; control surface positions ^aileron 
6 a and rudder 6 r ^; and time t. The full-scale range of the instruments is given in 
table HI. All the data were recorded on magnetic tape by an onboard recorder using a 
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pulse -code -modulation (PCM) recording system. Additional information on the data acqui- 
sition system can be found in reference 6. 

Data Preparation 

The data were initially recorded, digitized, and calibrated at Dryden Flight Research 
Center. A digital tape with the data in engineering units was sent to Langley Research 
Center. The acceleration data were corrected for instrument location. The difference in 
total pressure and static pressure measurements was assumed to be the dynamic pressure. 
Density was determined from the standard atmosphere tables for the measured pressure 
altitude and the airspeed was calculated from dynamic pressure. The indicated angles of 
attack and sideslip were corrected for the effects of aircraft angular rates. The linear 
velocities along the vehicle body axes were calculated from the airspeed, angle of attack, 
and angle of sideslip. All data were recorded on tape at the rate of 20 points per second. 
The tapes were then ready for use in the extraction program. 

RESULTS AND DISCUSSION 


Data for the flight conditions listed in table II were used with the mathematical model 
given in order to determine iteratively a set of aerodynamic derivatives for each of the 
flight conditions. The measured and computed time histories for each flight condition are 
shown in figures 2 to 6, with the measured data represented by dotted lines. The computed 
time histories shown are those attained after the differences between the measured and cal- 
culated time histories became small. The difference was considered small whenever 


R 


R 


i+1 


l R l 


< 0. 1, where R 


is the estimate of the error covariance matrix, as defined 


in the appendix, and i refers to the ith iteration. This inequality will be referred to as 
the cost function criterion. The figures show that in all tests, the computed time histories 
were in close agreement with the flight-record time histories. It should also be noted that 
the parameter values determined were consistent from run to run. Table IV gives the stan- 
dard deviations of the computed states from the measured states at convergence. The 
inverse square of each quantity in table IV was used as a diagonal term in the weighting 
matrix to obtain the fit of computed data to flight data shown in figures 2 to 6. The stan- 
dard deviation of each fit can be seen to be less than 3 percent of the full-scale instrument 
range, which was assumed to be the uncertainty in the measured data (see table III). For 
some quantities, the standard deviation was less than 1 percent of the full-scale measured 
quantity. 


The derivatives extracted for each flight condition (the derivatives which resulted Jin 
the computed time histories of figs. 2 to 6) are listed in table V along with their estimated 
standard deviations (Cramer-Rao lower bound). It should be noted that for M = 0.81 
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and = 8.2°, a single set of numbers is given which fits the data for both aileron input 
and rudder input. The single set was determined by using the values from the two runs 
that were least affected by correlation and that had the smallest standard deviations. 

As can be seen, this single set gave a good fit to all the data taken at M = 0.81 with 
a t = 8.2°. If the estimated standard deviations of the parameters were less than 10 per- 
cent of the values extracted, confidence in the values obtained was indicated. The param- 
eters which were not considered well determined were Cy , CY r , Cy fi , and C n g . 

The derivatives which were not considered well determined do not usually have a strong 
effect on the motion of the airframe. Therefore, it is difficult to excite the airframe so 
that there will be sufficient information in the data for a confident extraction of these 
derivatives. 


When analyzing the extraction results, the effects of correlation must also be con- 
sidered. If correlations between parameters are high, erroneous parameter values can 


result, and the inversion of the matrix 


l m M t 

i= 1 




(defined in the appendix) 


can be affected. The most obvious effects of correlation are two parameters assuming 
erroneous values but the fit to the data remaining good, or the cost function oscillating or 


diverging because the matrix 



is nearly singular (refs. 2 and 8). 


The correlation matrices for the aerodynamic parameters extracted are shown in 
tables VI(a) to VI(e). Although some correlations were high, none of the extracted parame- 
ter values seemed unreasonable, and convergence problems were encountered in only one 
run. 


For the run in which convergence problems occurred, a procedure discussed in ref- 
erence 2 was initially used to determine parameter values. This procedure first examined 
the covariance matrix to determine which parameters were potentially correlated. The 
covariance matrix for the run with = 13° and M = 0.66 indicated that C was 


correlated with C 7 and C; and that C n was correlated with C n and C n . For 

Lq ‘-n n fi_ 11 fl “n 


u p 


"P 


p 


this particular run, the correlation affected the convergence of the algorithm. To improve 
the convergence, either and C ng or C^, C^, C n ^, and c n were alter- 
nately held fixed. Initially and C ng were held fixed, and , C^, C n ^, and 

C n were allowed to change for several iterations. Then Ci , C 7 , C n , and C n 
n p & 6 /3 L p n /3 n p 
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were allowed to vary. This procedure was repeated 


were fixed, and C 7 - and C n 

°a °r 

until the change in the cost function was less than the criterion discussed earlier. 

Once a set of numbers was determined to have given a good fit to the data, all the 
parameters were made active. The fit showed almost no improvement, the parameter 
values changed less than 5 percent, and the system no longer diverged. Apparently the 
correlations were not strong enough to affect convergence when the initial guess at the 
parameter values was close to the converged values. This conclusion is supported by the 
fact that for the other runs no serious convergence problems were encountered, even though 
there were correlations for those runs which were as high as those for the rim in which 
convergence problems were initially encountered. The results shown in table VI(b) and in 
figure 3 are for the final part of the run when all the parameters were active. 

For comparison, values of some aerodynamic coefficients obtained from reference 6 
are shown in table VII. The numbers from reference 6 are for an altitude of 12. 19 km 
(40 000 ft), were transformed from stability to body axes, and were converted to a center 
of gravity located at 0.29c. Several of the derivatives that strongly affect the motion of 
the airframe were calculated using the method of reference 9, and these are also shown in 
table VII. 

Since there were some significant differences between the extracted values and the 
values from references 6 and 9, an effort was made to determine how well mathematical 
models using the various sets of parameters represented the airplane. The numbers from 
reference 6 that are given in table VII were put into the equations of motion, and for the 
proper Mach number and angle of attack, calculated responses are shown in figures 7 to 9. 
The fits to the flight data, although not as good as those obtained with the values from 
table V, are still reasonably good. 

In another effort to compare the different sets of parameter values, several of the 
sets were used to determine the stability characteristics of the mathematical model. These 
sets are designated A to F, and some of the conditions under which they were generated are 
given in table VIII. The comparison was made by use of an analysis from reference 10. 

The reader is cautioned that the approximate analysis given in reference 10 will not work 
for these cases unless all nonzero trim states are included in the equation of motion when 
the stability quartic is derived. 

The following equations were used in determining the damping ratios and periods of 
the Dutch roll: 
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Y -direction: 
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<p = p + r tan cos 0^. 


The damping ratios and the periods of the Dutch roll for various mathematical models are 
given in table Vin. Also given in table VIH are the periods and damping ratios of the Dutch 
roll mode calculated from the actual flight data (where there was sufficient free oscillation 
to make the calculations). As expected, when the period and damping ratios were calcu- 
lated using a mathematical model based on parameter values from reference 6, the damping 
ratios were less than when the extracted model was used. This conclusion is illustrated by 
figures 7 to 9. The analysis of the Dutch roll has demonstrated that several sets of param- 
eter values will give a fair fit to a set of flight data and similar stability characteristics; 
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however, for a given set of data, the extracted mathematical model gives the best fit to 
the data. 


CONCLUDING REMARKS 

Test data obtained during flights of the F-8C airplane at moderate to high angles of 
attack have been used to determine the lateral aerodynamic parameters of the airplane at 
four flight conditions. The tests were conducted with the airplane trimmed in a steady 
turn, with angles of attack approximately 9° and 13°, and with Mach numbers of 0.7 and 0.8. 

The extracted parameters were consistent from run to run and resulted in a fit to the 
flight data that varied by less than 3 percent of the full instrument range. The parameter 
values obtained were in fair agreement with values obtained from wind-tunnel and analytical 
methods. The adequacy of the extracted set of parameter values was further substantiated 
by showing that other sets of parameter values did not fit the flight data as well as the 
extracted set. However, period and damping ratios of the Dutch roll modes that were gen- 
erated by the parameter sets used for comparison were close to those generated by the 
extracted parameter set. 

Analytical and wind-tunnel studies have shown that several of the lateral aerodynamic 
parameters can vary with angle of attack over the angle-of-attack range tested. However, 
for each individual flight, the angle-of-attack variation from trim was so small that the 
linear aerodynamic model used was adequate to describe the motion of the airplane. There- 
fore, the existing parameter extraction program could be successfully used. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
December 1, 1976 
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APPENDIX 


PARAMETER ESTIMATION PROCEDURE 

The parameter estimation procedure used in this study is an iterative procedure, 
developed and discussed in reference 5. The procedure maximizes the conditional likeli- 
hood function L (aerodynamic parameters, weights, and initial conditions): 


L = 


1 

(2tt) 1 / 2 | Rl 1 / 2 


exp 




x 



where R is the estimate of the error covariance matrix, x is the vector describing the 
state of the airplane, T denotes the transpose, and -1 denotes the inverse. The states 
used in the likelihood function were v, p, r, 0, and ay. The calculated states c 
were determined by using the equations of motion previously introduced. In these equa- 
tions the longitudinal quantities u, w, q, and 9 were input directly into the equations 
from the flight data tape. The weighting matrix is R"*, where 


R = Diagonal 


N 1 

M*M-v)Krv) 


i= 1 


The result of maximizing the likelihood function is a parameter -updated equation which is 
given by 



r 

-1 

- 




N 

AC = 

N T 

i M (‘i) R -’ M ( t 1 ) 


1 M ('i) R '\ht - he) 


_i= 1 


_i= 1 J 


where C is the vector of aerodynamic coefficients to be determined, M is the matrix 
of sensitivities of the calculated states with respect to the unknown parameters. (See 


ref. 5.) The estimated parameter covariance matrix is 


J M(ti) T R" 1 M(ti) 


i=l 


The 


square root of each diagonal element of the covariance matrix (estimated standard devia- 
tion) is directly related to the uncertainty in the extracted parameters, and the off-diagonal 


13 



APPENDIX 


terms are used to indicate correlations between parameters. The program, whose devel- 
opment is discussed in reference 5, calculates the matrices and vectors required to gener- 
ate AC. The program then uses AC to change the aerodynamic coefficients iteratively 
until a fit to a set of flight data is obtained. The steps in the program operation are: 

(1) Choose values for the parameters to be identified. 

(2) Integrate the equations of motion using the current values of the aerodynamic 
parameters chosen, and get time histories of the states. 

(3) Compute the state covariance matrix R and the weighting matrix R"l. 

(4) Calculate the cost function, which is the determinant of R. 

(5) Integrate the set of differential equations for the sensitivities and then form the 
matrix M. 

(6) Form the maximum likelihood estimation equations for the parameter update 

AC. 

(7) Determine the new parameters using C = C curren t + AC. 

(8) Use the new C as the current value for the next iteration, and continue the pro- 
cess at step (2) until the cost function criterion is satisfied. 
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TABLE I.- GEOMETRIC CHARACTERISTICS OF THE F-8C AIRPLANE 


Fuselage length, m (ft) 16.52 (54.17) 

Wing: 

Area, m 2 (ft 2 ) 34.88 (375) 

Aspect ratio 3.4 

Span, m (ft) 10.88 (35.67) 

Mean geometric chord, m (ft) 3.59 (11.78) 

Vertical tail: 

Area, m 2 (ft 2 ) 10.14 (109) 

Aspect ratio 1.5 

Span, m (ft) 3.89 (12.75) 

Rudder area, m 2 (ft 2 ) 1.17 (12.56) 

Horizontal tail: 

Area, m 2 (ft 2 ) 8.68 (93.4) 

Aspect ratio 3.5 

Span, m (ft) 5.52 (18.1) 

Length (center of gravity to one-fourth of tail mean geometric 

chord), m (ft) 5.06 (16.6) 
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TABLE E.- FLIGHT CONDITIONS 


Run 

Nominal altitude 



Nominal 

Mach 

number 

Nominal 

deg 

Trim elevator 
deflection, 

5 e,t> 

deg 

Trim bank 
angle, 

h’ 

deg 

m 

ft 

-- 

1 

1 

34 000 

0.71 

9.2 

-8.13 

50 

2 

1 

34 000 

.66 


-9.74 


3 and 4 


KliljljB 

.81 

8.2 

-7.45 

-63 

5 

10 370 

1 

.78 

12.0 

-10.3 

70 



Moments of inertia 2. 

ivxac 

D 

h 

I 

Y 

1 

Z 

r xz 

kg 

slugs 

kg-m 2 

slug-ft 2 

kg- m2 

slug-ft 2 

kg-m 2 

slug-ft 2 

kg-m 2 

slug-ft 2 

9574.41 

656.06 

12 500 

9200 

118 000 

86 800 

124 000 

91 600 

4030 

2970 


a Any errors in assuming nominal conditions were no greater than 3 percent of the system uncertainty 
on the estimated values of the mass and inertias. 








TABLE m.- INSTRUMENT RANGES 3, 



Si 

v was calculated from v = V sin /3. Individual sensors are bas- 
ically more accurate than 3 percent of full scale; however, because of 
unknown errors, the effects of incompatibilities between measured states 
and processing errors, the system accuracy was assumed to be 3 percent 
of the full-scale range of the instrument for the data used during the 
extraction procedure. 
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TABLE IV.- STANDARD DEVIATIONS OF COMPUTED STATES 
FROM MEASURED STATES AT CONVERGENCE - 


a t’ 

deg 

M 

\ 

Standard deviation of - 




P, 

rad/ sec 

r, 

rad/ sec 

<P, 

rad 

a Y , 

g units 

m/ sec 

ft/sec 

9.2 

0.71 

0.3048 

1.0 

0.013 

0.0041 

0.068 

0.008 

13.0 

.66 

.4023 

1.32 

.018 

.0031 

.031 

.013 

a 8.2 

.81 

.2804 

.92 

.016 

.0026 

.019 

.010 

b 8.2 

.81 

.3200 

1.05 

.028 

.0022 

.036 

.019 

a 12.0 

.78 

.4755 

1.56 

.030 

.0040 

.048 

.011 


a Aileron input. 
^Rudder input. 
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TABLE V.- EXTRACTED VALUES OF PARAMETERS AND 
ESTIMATED STANDARD DEVIATIONS 



Extracted values (and standard deviations) for - 

Parameters 

«t = 9 - 2 ° 

«t = 13 ° 

Q 

ri- 

ll 

00 

iso 

0 

a t = 12° 


M = 0.71 

M = 0.66 

M = 0.81 

M = 0.78 






(a) 

(b) 

C Y e 

-0.80 

-0.80 

-0.80 

-0.74 

P 

(0.03) 

(0.04) 

(0.06) 

(0.035) 

CY D 

0.18 

0.17 

0.37 

0.13 

P 

(0.09) 

( 0 . 11 ) 

(0.12) 

(-0.15) 

c Y r 

0.45 

0.45 

0.45 

0.45 


(0.45) 

(0.7) 

(0.5) 

(-0.6) 

Cy 6 

0.11 

0.06 

0.06 

0.06 

°r 

(0.036) 

(0.03) 

(0.026) 

(Fixed) 

C» 

-0.12 

-0.112 

-0.14 

-0.11 


(0.003) 

(0.0013) 

(0.0015) 

(0.002) 

Cj 

-0.31 

-0.28 

-0.36 

-0.22 

L p 

(0.009) 

(0.0034) 

(0.004) 

(0.008) 

C lv 

0.51 

0.33 

0.37 

0.24 

c^r 

(0.027) 

(0.027) 

(0.03) 

(0.02) 

C Zg 

0.032 

0.028 

0.037 

0.03 

°r 

(0.0024) 

(0.0009) 

(0.00067) 

(Fixed) 

% 

-0.051 

-0.046 

-0.055 

-0.044 

6 a 

(0.0012) 

(0.0007) 

(0.00052) 

(0.0008) 

C n« 

0.102 

0.103 

0.11 

0.095 | 

p 

(0.003) 

(0.0018) 

(0.002) 

(0.003) 

C n 

-0.045 

-0.036 

-0.06 

-0.04 

P 

(0.001) 

(0.0045) 

(0.003) 

(0.01) 

C n 

-0.40 

-0.62 

-0.33 

-0.33 

1 r 

(0.026) 

(0.027) 

(0.02) 

(0.02) 

c n fi 

-0.12 

-0.13 

-0.115 

-0.115 

6 r 

(0.0027) 

(0.0013) 

(0.0007) 

(Fixed) 

C n fi 

-0.012 

-0.01 

-0.0084 

-0.0004 

6 a 

(0.0013) 

(0.0009) 

(0.00079) 

(0.001) 


a The values given were determined by taking a weighted average of the results 


from runs using aileron inputs only or rudder inputs only. 

b The vehicle was excited by using an aileron input only. 



TABLE VI.- CORRELATION MATRICES FOR EXTRACTED AERODYNAMIC PARAMETERS 


(a) a)- = 9.2° and M = 0.71, with both aileron and rudder inputs 



C Y 

Y /3 

Cy 

P 

Cy 

r 

Cy 6 

u r 

c i 

L (3 

c, 

L p 

c i 

L r 

C *6 

u r 

c? 

% 

S 

c ». 

S 

°r 


% 

1 

0.002 

-0.093 

-0.2 

0.088 

-0.035 

0.072 

0.095 

0.064 

0.13 

0.15 

-0.19 

-0.19 

0.19 

% 

0.002 

1 

-0.58 

-0.13 

-0.034 

-0.01 

0.015 

-0.024 

-0.064 

-0.060 

-0.066 

-0.045 

0.063 

-0.021 

Cy 

i Y 

-0.093 

-0.58 

1 

0.3 

0.06 

0.03 

0.18 

0.003 

0.046 

-0.032 

-0.044 

-0.16 

-0.04 

-0.06 

% 

-0.2 

-0.13 

0.3 

1 

0.016 

0.006 

0.014 

0.021 

0.025 

-0.038 

0.009 

-0.011 

-0.002 

-0.02 

% 

0.088 

-0.034 

0.06 

0.016 

1 

0.63 

-0.37 

-0.46 

0.70 

0.73 

0.51 

-0.61 

-0.59 

0.59 ' 

% 

-0.035 

-0.01 

0.03 

0.006 

0.63 

1 

-0.54 

-0.35 

0.89 

0.40 

0.45 

-0.36 

-0.13 

0.44 

Cl 

, L r 

0.072 

0.015 

0.18 

0.014 

-0.37 

-0.54 

1 

0.50 

-0.45 

-0.43 

-0.35 

-0.19 

0.13 

-0.41 


0.095 

-0.024 

0.003 

0.021 

-0.46 

-0.35 

0.50 

1 

-0.35 

-0.30 

0.018 

0.025 

0.16 

-0.10 


0.064 

-0.064 

0.046 

0.025 

0.70 

0.89 

-0.45 

-0/35 

1 

0.46 

0.48 

-0.44 

-0.22 

0.40 

C «i3 

0.13 

-0.060 

-0.032 

-0.038 

0.73 

0.40 

-0.43 

-0.30 

0.46 

1 

0.76 

-0.63 

-0.76 

0.86 

° n P 

0.15 

-0.066 

-0.044 

0.009 

0.51 

0.45 

-0.35 

0.018 

0.48 

0.76 

1 

-0.68 

-0.53 

0.93 

c n r 

-0.19 

-0.045 

-0.16 

-0.011 

-0.61 

-0.36 

-0.19 

0.025 

-0.44 

-0.63 

-0.68 

1 

0.68 

-0.64 

S r 

-0.19 

0.063 

-0.04 

-0.002 

-0.59 

-0.13 

0.13 

0.16 

-0.22 

-0.76 

-0.53 

0.68 

1 

-0.61 

S a 

0.19 

-0.021 

-0.06 

1 -0.02 

0.59 

0.44 

-0.41 

-0.10 

0.40 

0.86 

0.93 

-0.64 

-0.61 

1 




























































TABLE VH.- PARAMETER VALUES FROM OTHER SOURCES 


Parameters 

Values from ref. 6 

Values obtained 

Case E 

Case F 

using ref. 9 

8 ° 

a t = 12 ° 

a t = 8.2° 

M = 0.8 

M = 0.8 

M = 0.81 

Cy 

Y /3 

-1.00 

-1.00 

-0.82 

Cy 

P 

0.175 

0.18 


Cy 

r 

0.45 

0.45 


Cy 

Yg r 

0.39 

0.195 


% 

-0.126 

-0.092 

-0.115 

% 

-0.240 

-0.085 

-0.336 

% 

0.057 

0.070 


c \ 1 

0.014 

0.007 


C *5a 

-0.060 

-0.Q60 


C n 

n /3 

0.143 

0.092 

0.17 

c n 

P 

-0.032 

-0.035 


Cn r 

-0.30 

-0.33 

-0.17 

Cn 6 r 

-0.104 

-0.052 


c n« Q 

-0.004 

-0.0005 




TABLE Vm.- PERIODS AND DAMPING RATIOS FOR THE DUTCH ROLL MODE 


Case 

M 

a t> 

deg 

P f> 

sec 

P c’ 

sec 

*f 

5c 

Conditions 

A 

0.71 

9.2 

2.2 

2.6 

0.2 

0.23 

Parameters extracted from flight data 
as given in table V. 

B 

.66 

13 

2.2 

2.5 

.2 

.23 

Parameters extracted from flight data 
as given in table V. 

C 

.81 

8.2 

2.2 

2.3 

.2 

.20 

Parameters extracted from flight data 
as given in table V. 

D 

.78 

12 


2.2 


.19 

Parameters extracted from flight data 
as given in table V. 

E 

a .8 

8 

2.2 

2.0 

.2 

.17 

Parameters from reference 6 as given 
in table vn. 

F 

a .8 

12 


2.2 


.12 

Parameters from reference 6 as given 
in table VII. 


a From reference 6 model. 
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16.52 m (54.17 ft) 













t, sec 


Figure 3.- Comparison of measured data with time histories computed by 
using parameters given in table V for flight data taken at M = 0.66 
and at = 13°, with both aileron and rudder inputs. 
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Figure 7.- Comparison of measured data of figure 4 with the time histories computed 
by using the parameters of reference 6 as given in table VII for Case E. 
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